arXiv: 1503.05102vl [gr-qc] 17 Mar 2015 


Homoclinic accretion solutions in the Schwarzschild—anti-de Sitter spacetime 
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The aim of this paper is to clarify the distinction between homoclinic and standard (global) Bondi- 
type accretion solutions in the Schwarzschild-anti-de Sitter spacetime. The homoclinic solutions 
have recently been discovered numerically for polytropic equations of state. Here I show that they 
exist also for certain isothermal (linear) equations of state, and an analytic solution of this type 
is obtained. It is argued that the existence of such solutions is generic, although for sufficiently 
relativistic matter models (photon gas, ultra-hard equation of state) there exist global solutions that 
can be continued to infinity, similarly to standard Michel’s solutions in the Schwarzschild spacetime. 

In contrast to that global solutions should not exist for matter models with a non-vanishing rest- 
mass component, and this is demonstrated for polytropes. For homoclinic isothermal solutions I 
derive an upper bound on the mass of the black hole for which stationary transonic accretion is 
allowed. 

PACS numbers: 04.20.-q, 04.40.Nr, 98.35.Mp 


In a series of 3 papers m Karkowski, Malec, and I 
have discussed stationary, Bondi-type accretion flows in 
Schwarzschild-de Sitter and Schwarzschild-anti-de Sit¬ 
ter spacetimes. We derived analytic solutions for three 
“isothermal” equations of state of the form p = fee, where 
p is the pressure, and e denotes the energy density, with 
k = 1/3,1/2, and 1. For polytropic equations of state 
of the form p = Kp^, where p is the baryonic density 
and K and F are constant, we have computed numerical 
solutions. 

While in the sector of the positive cosmological con¬ 
stant A all these solutions behave in a similar fashion, 
their behavior differ significantly for negative A. In par¬ 
ticular, isothermal transonic solutions with k = 1/3,1/2, 
and 1 cover the entire space outside the black hole 
horizon, i.e., they can be continued to arbitrarily large 
radii (similarly to standard Bondi-type solutions in the 
Schwarzschild case mm- The situation is different for 
polytropic solutions, and this behavior is illustrated in 
Figs. 2 and 3 of [2]. When a polytropic solution is con¬ 
tinued outwards, one encounters a finite radius r, where 
the derivative of the radial velocity component dru'" di¬ 
verges. 

This behavior corresponds to the existence of homo¬ 
clinic orbits in the phase diagram of the radius r versus 
the radial component of the velocity u”. Given the equa¬ 
tions describing stationary Bondi-type accretion one can 
construct a dynamical system with the phase portrait 
consisting of the graphs of different solutions M”(r) de¬ 
scribing the accretion flow [5]. The homoclinic orbit of 
this dynamical system consists of two transonic solutions 
— the solutions that pass through a saddle-type critical 
(fixed) point, at which the local value of the speed of 
sound is equal to the modulus of the three velocity of the 
fluid (the so-called sonic point). 

Given the above examples, one could have an impres¬ 
sion that those homoclinic solutions are somehow essen¬ 
tially related to the assumed polytropic equation of state 
(as opposed to the isothermal form). This is not quite 
true, and the aim of this paper is to clarify these issues. 


It is possible to show, and I do that in Sec. |TT] of this 
paper, that for isothermal equations of state a power-law 
asymptotic behavior is admitted only if /c > 1/3. This 
suggests the existence of homoclinic solutions also for 
isothermal equations of state with k < 1/3, and indeed 
an analytic example of such a solution can be given for 
k = 1/4 (Sec. [ml). 

On the other hand, one can show (Sec. that global 
(extending to infinity) solutions do not exist for poly¬ 
tropic equations of state, irrespective of the assumed 
values of the speed of sound. Interestingly, such solu¬ 
tions are not permitted due to the term in the expression 
for the specific enthalpy that represents the contribution 
from the rest-mass of gas particles. 

The above observations suggest an analogy with the 
properties of radial time-like and null geodesics in the 
Schwarzschild-anti-de Sitter spacetime [7]. For time-like 
geodesics the term proportional to in the effective po¬ 
tential forbids the motion of a particle with a finite energy 
to infinity. In contrast to that, a massless particle can 
travel to arbitrary large radii. 

It seems quite plausible that the occurrence of homo¬ 
clinic solutions is generic, and it is characteristic for non- 
relativistic matter, where the energy density has a non¬ 
vanishing contribution from the rest-mass of the gas par¬ 
ticles. 

Somewhat on the margin of the above considerations I 
show that for isothermal equations of state with fc < 1/3 
transonic accretion solutions exist only for sufficiently 
small black holes. The appropriate limit on the mass 
of the black hole is derived in Sec. IIVI 

Relativistic Bondi-type accretion flows were inves¬ 
tigated in many recent papers. For instance, self- 
gravitating flows were analyzed in [HIS], and also in m- 
Radiation transfer was included in [111 I12j . Among very 
recent developments one should also notice [T3]. There 
are also papers dealing with Schwarzschild-anti-de Sitter 
spacetimes. In [14] Amani and Farahani discussed phan¬ 
tom accretion onto Schwarzschild-anti-de Sitter black 
holes, basing on the idea presented in m- Accretion 
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onto Schwarzschild-anti-de Sitter black holes is also con¬ 
sidered in [16]. 


I. NOTATION AND EQUATIONS 


In this paper I assume the gravitational system of units 
with c = G = 1, and the signature of the metric tensor 

+j +: +)• 

The Schwarzschild-anti-de Sitter metric in standard 
polar coordinates (t, r, 0, (p) has the form 


ds^ 



-br^ (d0^ -b 



sin^ , 


dr'^ 


( 1 ) 


where the cosmological constant A is negative, and m de¬ 
notes the mass of the black hole. Alternatively, one can 
work in suitably chosen Eddington-Finkelstein-type co¬ 
ordinates, that are regular at the horizon. An interested 
reader is referred to |5] . 

The horizon of the black hole is located at 
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- ar sinh 


(^3to\/|A| 


(a real positive root of the equation 1 — 2m/r — Ar^/3 = 

0 ). 

The motion of the fluid is described by the conservation 
laws 


V^T'^"=0, V^(pu'^)=0, (2) 


both cases p, e, h, and p are related by simple formulas, 
provided that the flow of the gas is smooth, i.e., there 
are no shock waves or contact discontinuities. For the 
isothermal equations of state one has 

e = Glpl+^ h = {l + k)Cip\ (4) 

where Gi is a constant. In the case of the polytropic 
equation of state, it is easy to show that 

e = Gap -b h = C 2 + (5) 

The choice of the constant G 2 is subject to the physi¬ 
cal interpretation. In the following I assume standard 
normalization with G 2 = 1. By setting G 2 = 0 one can 
recover Fqs. 

It is also important to stress that the conservation law 
= 0 in Fqs. ([^ is, de facto, the definition of p. 
This quantity can have different physical interpretation, 
depending on the context. For the gas of a conserved 
number of massive particles, the density p can be ex¬ 
pressed as the number density of the particles times the 
mean rest-mass of the particle. On the other hand, the 
dynamics of the photon gas with the equation of state 
p = e/3 is described by equations ^= 0 only. In 
this case it is also possible to introduce a function p sat¬ 
isfying the conservation law = 0, provided that 

the flow of the gas is smooth. Such a function has a nat¬ 
ural interpretation of the specific entropy (note that it is 
not conserved across possible discontinuities in the flow). 

By differentiating Fqs. (|^ with respect to r one can 
obtain the following expression for du^/dr: 


where 


= {e + p)uf^u''+pg^^’' 


du^ 

dr 


2u- c2 [1 - ^ - A+ (uO"] - f + 

r (u’-)2 _ c2 [I - 2 ^ - Ar2 + [u^f] 


is the energy-momentum tensor. Here e, p, and p de¬ 
note the energy density, the pressure, and the baryonic 
density, respectively. In what follows we introduce the 
specific enthalpy h = {e + p)/p (note that the specific 
enthalpy is sometimes defined as h = {e+p)/p—\, which 
is different form the convention used in this paper). The 
symbol denotes the components of the metric tensor; 
are the components of the four-velocity of the fluid. 

For spherically symmetric and stationary flows = 
= 0, and all quantities appearing in Fqs. ([^ are func¬ 
tions of r only. Assuming that the solution is smooth, 
one can integrate Fqs. ([^. This yields 


h\ 1 -r2 -b (u’’)2 = const, r'^pu'' = const. 

V r 3 

( 3 ) 

The above equations constitute a starting point for the 
discussion of the Bondi-type accretion solutions. 

In this work I deal with two classes of equations of 
state. The isothermal equations of state p = fee, where 
0 < /c < I is a constant HZ], and polytropic equations of 
state p = Kp^, where K > 0 and T > 1 are constant. In 


where dh/dp = he/./p, and is the speed of sound. One 
has c/ = k for isothermal equations of state, and = 
(T —1)(1 — 1 // 1 ) for polytropes. Formally, this constitutes 
a “step backwards” in the process of finding of solutions, 
but it gives an insight into their structure. Let I — l{r) 
be a function of r such that 


dr f t 2 

= fi{r, u^)- 


1 - 


2m A o 

- r^ 

r 3 




( 6 ) 


Then 




^/2(r,u”). 


A 
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+ {vdf 


m 



( 7 ) 


Equations Q and 0 can be treated as a dynamical sys¬ 
tem, whose phase portrait consists of the graphs of u” 
versus r, or more precisely, the graphs of u”(r) belong to 
the orbits of Fqs. § and 0- ^ ^ 

The dynamical system defined by Fqs. (|^ and 0 
has critical (fixed) points (r*,u”) where /i(r*,u”) = 
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/ 2 (r *,0 = 0, that is 
{<f - cL 


2?71 iV 2 / r \ 2 

1- o + {<? 

r* c> 


9 2??7- iV 9 / T \ 9 

"" 1- 

r* 3 


m 


A 2 

+ -T^rl 


0 , ( 8 ) 
0. (9) 


Here, and in what follows, the quantities referring to the 
critical point will be denoted with an asterisk. 

In order to find critical points one has to specify the 
equation of state. It is, however, clear that 


2 _ A 2 

2r* e''*’ 


( 10 ) 


and also 

m _ A 2 
„2 2r, 6 ' * 

1 3m A 2’ 

2r. 2 ' * 

Thus, one can easily show that at the critical point 
cL = {9rrV''v'')t = {VrV'")^, 

where u’' = ju* is the radial component of the three- 
velocity. Thence it is common to use the term “sonic 
point” instead of “critical point”. For A = 0, i.e., in the 
Schwarzschild case, and for standard equations of state, 
(r*,u”) is a saddle point. Thus, there is an accretion 
solution passing through the sonic point. In the following 
I will reserve the term “sonic point” for a critical saddle 
point, as opposed to other types of critical points. 

The existence and number of critical points will be 
discussed in Sec.j^for the isothermal equations of state. 


II. ASYMPTOTIC BEHAVIOR OF 
ISOTHERMAL SOLUTIONS 

It is quite illuminating to start with an anaNsis of 
the asymptotic behavior of isothermal solution^ For 
isothermal equations of state Eqs. ^ yield 

l_2m_A2^ 2^^^4fe|^.|2fc^ (11) 

r 3 

where A is a positive constant. Let us assume an asymp¬ 
totic expansion of the form u” ~ Hr“. One gets 

= AB2fcr2/c(2-sa)^ 

r 3 

There are two ways in which this equation can be satisfied 
in the leading order for r —> oo. In both cases the term 
on the right-hand side has to cancel with the leading 
term at the left-hand side. If the leading order term on 


^ A part of the analysis presented in this chapter was pointed to 
me by Jerzy Knopik. 



FIG. 1. Solutions obtained for the equation of state p = e/4 
and sample metric parameters m = 1, A = —5/100. The solid 
line depicts transonic solutions. The vertical line denotes the 
position of the horizon of the black hole. The graph shows 
solutions for A' = AjA^ = 0.98, 1, and 1.02. 


the left-hand side is then 2k(2 + a) = 2, and 

2 > 2a (the condition that the term is not the 

leading one). If the leading order term is we have 

2k{2 -I- a) = 2a, and 2a > 2. These two possibilities lead 
to a = (1 — 2k)/k and a = 2A:/(1 — fc), respectively. In 
either case k > 1/3. For fc = 1/3 the two asymptotics 
coincide; one has a = 1. It can be seen from the examples 
given in [2] that the two exponents a correspond to the 
branches of the solution that are asymptotically subsonic 
and supersonic respectively (note that for 1/3 < fc < 
1 we have (1 — 2k)/k < 2fc/(l — k)). If fc = 1/3, the 
two branches still exist, but with the same asymptotic 
behavior. 

It is clear that for fc < 1/3 asymptotic behavior of 
the form u” ~ Br°‘ is not admitted. This suggests that 
for fc < 1/3 global solutions do not exist. The following 
example confirms this expectation. 


HI. SOLUTION FOR p = e/4 

An explicit homoclinic solution can be found for the 
equation of state p = e/4. For this equation of state 
Eq. can be written as 

f(12) 

Here I will be mainly interested in transonic solutions. In 
this case the value of A should be determined by requiring 
that the solution passes through the sonic point (a saddle 
critical point of the dynamical system defined by Eqs. ([^ 
and 0). The location of the sonic point r* is given by 

^ i /^l _ ^ _ (^^2\ 

2r* 6 * 2r* 2 
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or 


Arl + 6r* — 21m = 0 


(cf. Sec. IV). The above polynomial equation has real 


positive solutions only for m^ < 32/(441|A|), i.e., when 
its discriminant is not positive. This follows from a very 
simple reasoning; it is presented in Sec. |IV| without re¬ 
stricting to A: = 1/4 only. For < 32/(441|A|) the 
location of the sonic point is given by 


2^2 

r* = ,— cos 


^ I ^arcco'’ 

- + -arccos|^ 4 ^ 


The square of the radial component of the velocity at the 


sonic point is given by Eq. (101; the value of the constant 
A appearing in Eq. (12) that corresponds to transonic 


solutions can be expressed as 


A — A^ = — 
r* 


m 

2r* 


A . 

- r‘ 

6 ^ 


(13) 


Equation (121 is a quartic polynomial equation in 


and it can be solved exactly. Let 


^ 4 j 33 


1 - 


-r 
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by the second equation in Eqs. (3). The expressions for 
e and h are obtained from Eqs. (4|. 

Examples of solutions for the equation of state p = e/4 
are depicted in Eig. 1. Here m = 1, A = —5/100, and the 
solutions are plotted for A' = A/A^ = 0.98,1, and 1.02, 
respectively. 


IV. CRITICAL POINTS IN THE ISOTHERMAL 
CASE 

The discussion of the asymptotic behavior of isother¬ 
mal solutions given in Sec. [IT] and the solutions shown 
in Fig. 1 suggest that for isothermal equations of state 
p = ke with 0 < /c < 1/3 one should expect the existence 
of two critical points on the phase diagram of (r, m’'), at 
least for a certain range of parameters. These should be 
a saddle point belonging to the homoclinic orbit, and a 
center enclosed by this orbit. I will now argue that this 
is indeed the case. In addition, a simple reasoning gives 
a restriction on the mass of the black hole allowing for 
transonic accretion. 

For isothermal equations of state p = ke Eqs. (|^ and 
@ yield 


m 

2r* 


- = fc ( 1 - 


3m A . 

-C 

2r* 2 * 


where A = A* is given by Eq. (13). The domain of the which is equivalent to the cubic equation 


transonic solution is given by the condition A > 0. Let 
R denote the largest root of the equation A = 0, and let 
us set 


Y -y 


1± 


' Ar 

V2yi 


- 1 


y = 







where in the expression for y we choose real roots. Then 
the two branches of the transonic solution can be written 
as 

0<r<r., 

\x_, <r < R, 

(this branch is subsonic outside r*) and 
' ' \X+, r,<r<R, 

for the branch that is supersonic outside r*. For values 
of the constant A other than that given by Eq. (13) the 


expressions X± give remaining accretion solutions, that 
are not transonic. 

Once the expressions for rt’’ are found, all other quan¬ 
tities can be trivially computed. The density p is given 




rt -\- kr^ — [ fc -I- ^ ) = 0. 


Dividing by ^ — fc) one can write the above equation 


as 


where 


fir*) = rt + Zpr* -f 2q = 0, 


(14) 


2k 

A(l-3fc)’ 


3m(3k -I- 1) 
'2A(l-3/c) ■ 


For A < 0 and 0 < fc < 1/3 the above equation has 
a real negative root. This is a trivial consequence of the 
fact that /(O) = 2q > 0 and limf.._>._oo f{r*) = — 00 . 
Thus, a real and positive root can exist only if the dis¬ 
criminant W = +p^ is nonpositive. A straightforward 

calculation shows that this condition is equivalent to 


m^ < 


32fc3 


9(l-3fc)(l-f 3/c)2|A|’ 


(15) 


which can be interpreted as an upper bound on the black 
hole mass allowing for the transonic accretion. In this 
case Eq. (14) has three real roots, two of them being 
positive. This last statement follows directly from the 
analysis of the complex roots in the Cardano formula 


r* = \ - 


\J-q + i\/W\ + \J-q-i\/W\- 
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Conversely, for fc > 1 /3 and A < 0 the discriminant W 
is always positive, and there is just one real root. This 
root is positive: in this case /(O) < 0 and, of course, 
limr,_>oo f(r*) = + 00 . 

In principle, the analysis of the critical points can be 
pursued further by computing the Jacobians of the right- 
hand side of Eqs. ^ and Q, i.e.. 


^ % 
4 IS 


dr du^ 


(16) 


at critical points, and analyzing their eigenvalues. This 
can be easily done, say with Wolfram Mathematica. One 
can show that the critical points are either a saddle (Ja¬ 
cobian (161 has real eigenvalues with different signs) or 
a center (eigenvalues of (161 are imaginary). More pre¬ 
cisely, it can be shown that the eigenvalues are of the 
form ±\/Y, but the resulting expression for Y is lengthy, 
and the sign of Y is not immediately clear. A case-by- 
case study shows that the critical point with smaller r* is 
a saddle point [Y >0), and the critical point with larger 
r* is a center {Y < 0), as expected. For the data depicted 
in Fig. 1, i.e., k = 1/4, m = 1, and A = —5/100, the crit¬ 
ical points are located at r* = 4.05608 and r* = 8.34795. 
In the first case the eigenvalues of (16) are ±0.391511. 
In the second they are ±0.68955b 

The mass limit (15) is illustrated in Fig. 2. It shows 
a sequence of transonic solutions obtained for increasing 
values of the black hole mass and the equation of state 
p = e/4. The solid line in this graph depicts the limiting 


case with the mass m = rriL = -\/32/(441|A|). Note that 
the “homoclinic loop” gets smaller and smaller with the 
increasing mass, and it disappears for m = mz,. 


V. POLYTROPIC SOLUTIONS 


Sample numerical solutions with polytropic equations 
of state were given in [5], exhibiting characteristic ho¬ 
moclinic behavior. Here I note that no global polytropic 
solutions exist (i.e., Fqs. © cannot be satisfied asymp¬ 
totically), irrespective of the assumed parameters of the 
polytropic equation of state. This can be seen immedi¬ 
ately, if one attempts to repeat the reasoning of Sec. [IT] 
From Fqs. and ([^ one has 


h 



2m 

r 



TK 

r-1 


p 


r-i 



^2 ± ( u ”)2 = 

O 

^r2 ± (u’’)2 = const. 

O 


The square root in the above formula behaves asymptot¬ 
ically at least as 0{r), due to the A term. The only way 
in which this divergent behavior could be cancelled (to 
yield the constant appearing on the right-hand side) is to 
have h vanishing sufficiently fast. This is clearly impos¬ 
sible since h > 1. Of course, the unity in the expression 
for h comes from the p term in the formula for the energy 


ImJ 



FIG. 2. Transonic solutions obtained for the equation of state 
p = e/4 and A = —5/100. The graphs correspond to fonr 
values of the black hole mass: m = 0.85m_L, m — O.Omz, 
m = 0.95m_L, and m = mz, where mz = -y/32/(441|A|) is the 
maximal mass for which a sonic point exists. 


density e, i.e., a contribution due to the rest-mass of the 
particles of the gas. 


VI. DISCUSSION 


A very natural interpretation of the above results is to 
say that the existence of the global (asymptotic) solutions 
of the Bondi-type accretion in the Schwarzschild-anti-de 
Sitter spacetime is not directly related to the algebraic 
form of the equation of state. Global solutions exist for 
relativistic matter models that can be associated with 
the gas of massless particles. In this context it is im¬ 
portant to note that the equation of state p = a 
limiting case among isothermal equations of state, is a 
simple well-known model of the photon gas. This inter¬ 
pretation is also confirmed by the observation that the 
term that forbids asymptotic solutions in the case of poly¬ 
tropic equations of state is directly connected with the 
non-vanishing rest-mass of the gas particles. 

Of course, the simplest (and perhaps slightly naive) 
explanation would be to refer to the behavior of radial 
null and time-like geodesics in the Schwarzschild-anti- 
de Sitter spacetimes. They are analyzed in detail in [7]. 
Time-like radial geodesics are described by 



where r denotes the proper time, and E is constant (the 
energy). Here 1 — 2mfr — Ar‘^/3 plays a role of the ef¬ 
fective potential. Since it diverges asymptotically as r^, 
the range of the motion of the particle with finite energy 
E is limited. In terms of the coordinate time the above 
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equation can be written as 



The equation describing null geodesics has the form 



meaning that a photon can travel to infinity. 

It is also interesting to note that a similar homoclinic 
behavior was observed in the spherically symmetric mod¬ 
els of Bondi-type accretion on Reissner-Nordstrdm black 
holes m- In this case it is caused by the electric charge 


term, and the homoclinic loop is located inside the hori¬ 
zon of the black hole. 

As a by-product of the analysis presented in this note I 
obtained a restriction on the maximum mass of the black 
hole that allows for the transonic accretion of isothermal 
fluids with 0 < A: < 1/3. This bound is given by a re¬ 
markably simple formula, and it is discussed in Sec. EYl 
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